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THE MASSES, ANCESTORS AND DESCENDENTS OF EXTREMELY RED OBJECTS: 
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ABSTRACT 

Wide field near-infrared (IR) surveys have revealed a population of galaxies with very red optical— IR 
colors, which have been termed "Extremely Red Objects" (EROs). Modeling suggests that such red 
colors (R — K > 5) could be produced by galaxies at z ^ 1 with either very old stellar populations or 
very high dust extinction. Recently it has been discovered that EROs are strongly clustered. Are these 
objects the high-redshift progenitors of present day giant ellipticals (gEs)? Are they already massive at 
this epoch? Are they the descendents of the z ~ 3 Lyman Break Galaxies (LBG), which have also been 
identified as possible high redshift progenitors of giant ellipticals? We address these questions within 
the framework of the Cold Dark Matter paradigm using an analytic model that connects the number 
density and clustering or bias of an observed population with the halo occupation function (the number 
of observed galaxies per halo of a given mass). We find that EROs reside in massive dark matter halos, 



with average mass (M) > 10 n 100 m . 



n ioo ■ The occupation function that we derive for EROs is very similar 
to the one we derive for z — 0, L > L. t early type galaxies, whereas the occupation function for LBGs 
is skewed towards much smaller host halo masses ((M) « 10 n ~ 12 h^Q Mq). We then use the derived 
occupation function parameters to explore the possible evolutionary connections between these three 
populations. 

Subject headings: galaxies: high-redshift — galaxies: halos — galaxies: evolution — galaxies: 
fundamental parameters — theory: dark matter 



1. INTRODUCTION 

The first generation of deep near-IR surveys turned up 
a population of objects not represented in optical sur- 
veys (Elston, Rieke, & Rieke 1988; Elston, Rieke, & Rieke 
1989; McCarthy, Persson, & West 1992; Hu & Ridgway 
1994; Cowie et al. 1994; Moustakas et al. 1997; Barger 
et al. 1999; Thompson et al. 1999). These objects are re- 
markable for their very red optical— infrared colors (e.g. 
R — K > 6), and hence gained the moniker "Extremely 
Red Objects" (EROs). Using cither empirical templates 
or stellar population modeling, it is found that these ex- 
treme colors could be produced in galaxies at z ,> 1 with 
either a very old and quiescent stellar population or very 
large dust extinction (see e.g. Firth et al. 2002, Fig. 1). 
At these redshifts, the observed magnitudes K ^ 19 im- 
ply large luminosities, L > L*. Taken together, these clues 
hint that EROs could be the elusive high redshift progen- 
itors of present-day giant ellipticals, already massive and 
with evolved stellar populations at an epoch when the total 
age of the Universe is only 4 — 5 Gyr. If correct, this clearly 
has dramatic implications for theories of galaxy formation. 
In particular, it casts in sharp contrast the hierarchical, 
Cold Dark Matter (CDM) based scenario, in which ellip- 
ticals form in gas-rich mergers, and the (phenomenologi- 
cal) "monolithic collapse" scenario (Eggen, Lynden-Bcll, 
& Sandage 1962), in which ellipticals form at high redshift 



and evolve passively. 

For several years, further progress in the interpretation 
of these objects was hindered by the difficulty of obtain- 
ing reliable measurements of their number densities and 
redshifts. The surface densities of EROs measured in the 
early fields varied widely (e.g. Cowie et al. 1994; Mous- 
takas et al. 1997; Thompson et al. 1999), perhaps not 
surprisingly in view of the small volumes that had been 
probed. In addition, the faintness and apparent lack of 
dramatic spectral features has made a systematic spectro- 
scopic study of a large sample of EROs impractical. Heroic 
spectroscopic efforts directed at a few objects secured sev- 
eral high signal-to-noise-ratio spectra, demonstrating that 
some of the EROs indeed contain predominantly old stel- 
lar populations with very little recent star formation and 
are at redshifts of z ~ 1 (Spinrad et al. 1997; Dunlop 
et al. 1996). Some, however, appear to be heavily dust 
enshrouded starburst galaxies, similar to local ULIRGs, 
also at redshifts z ~ 1 (Graham & Dey 1996; Smail et al. 
1999). 

Recently, the completion of portions of several wide field 
surveys with multi-band, optical to near-IR photometry 
(e.g. Daddi et al. 2000; Firth et al. 2002; Chen et al. 2002; 
Olding et al. 2002; Martini 2001) has led to several in- 
teresting developments. The area surveyed is finally large 
enough that the global average surface density of EROs 
is more robustly measured (though still with some uncer- 
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tainty) to limiting magnitudes K « 19 — 20. These surveys 
have also shown that the EROs are strongly clustered on 
the sky, with an angular clustering amplitude ~ 10 times 
larger than the overall population at the same limiting 
magnitude. To obtain estimates of the real-space clus- 
tering, one needs redshift information. Obtaining large 
numbers of spectroscopic redshifts has proven to be ex- 
tremely difficult, but photometric redshifts indicate that 
EROs selected at, for example, I — H > 3, lie in a rela- 
tively narrow redshift range with a mean of z ~ 1.2 (Firth 
et al. 2002). This implies a real-space correlation length 
of r ~ 9 - lOhj"^ Mpc (McCarthy et al. 2001; Firth et al. 
2002; Daddi et al. 2001, precise values depending on the 
color and magnitude cut), comparable to or larger than 
z = luminous early type galaxies. 

Gradual progress is also being made in determining the 
break-down of the ERO population in terms of the two 
primary classes of "old" and "dusty" objects, using ju- 
diciously chosen broadband photometry (e.g. Pozzetti & 
Mannucci 2000), spectroscopy (Cimatti et al. 2002), and 
morphology (Moriondo et al. 2000). From all of these 
measures, the fraction of "old" objects is consistent with 
50±20%, although the contribution from dusty objects 
may be larger at faint magnitudes. Based on a small 
sample of 35 EROs with spectroscopic redshifts, Daddi 
et al. (2002) found that the "dusty" objects are signifi- 
cantly less clustered than the "old" objects. The correla- 
tion length for the old EROs is constrained to be in the 
range r ~ 5.5 — 16h^ Q Mpc, consistent with previous es- 
timates based on the full photometric samples. 

Knowledge about the clustering properties of these ob- 
jects provides important clues about their nature and their 
connection with other populations. In the Cold Dark Mat- 
ter (CDM) paradigm, massive halos correspond to rare 
peaks in the density field and are therefore strongly clus- 
tered. Thus we expect a correlation between halo mass 
and clustering amplitude, and the description of this cor- 
relation by analytic theories (Mo & White 1996; Shcth 
& Tormen 1999) agrees quite well with the results of N- 
body simulations. If the values of cosmological parameters 
are determined by other means, then in principle we can 
use the observed clustering properties of any population 
to constrain the possible masses of the host dark matter 
halos. In practice, an added complication is the unknown 
"occupation function" of the observed population; i.e. the 
number of observed galaxies per dark matter halo as a 
function of halo mass. For any number- or pair-weighted 
statistic like the correlation function, this effectively gives 
different weights to halos of different masses. The occu- 
pation function in essence provides the link between an 
observed population and the theoretically tractable popu- 
lation of dark matter halos. 

The occupation function is determined, in reality, by 
the host of interwoven astrophysical processes controlling 
gas cooling, star formation, feedback, and so on. It may 
be computed from first principles using detailed models 
of galaxy formation, i.e. semi-analytic models (Somcrvillc 
et al. 2001a; Benson et al. 2001) or hydrodynamic simula- 
tions (White et al. 2001). Numerous studies have run these 
"forward" models to fill dark matter halos with galaxies, 
and then evaluated whether the resulting number densities 
and clustering properties match observations for various 



populations (e.g. Wechsler et al. 2001; Somerville, Pri- 
mack, & Faber 2001b; Katz, Hernquist, & Weinberg 1999). 
Alternatively, an empirical approach may be taken, which 
bypasses the necessity of making any assumptions about 
the physics of galaxy formation. If a specific functional 
form for the occupation function is adopted, its parame- 
ters can be constrained by appealing directly to observa- 
tions. For example, Peacock & Smith (2000) and Marinoni 
& Hudson (2002) do this using the observed luminosity 
function of groups, and Bcrlind & Weinberg (2001) use 
the clustering properties of nearby galaxies. 

In this paper, we make use of a similar empirical ap- 
proach, using a simple analytic model for galaxy clus- 
tering, based on the formalism introduced in Wechsler 
et al. (2001) and Bullock, Wechsler, & Somerville (2001, 
BWS02). In this model, we parametrize the occupation 
function using a simple functional form, a power law char- 
acterized by three parameters: a minimum mass, a slope, 
and a normalization. This form may be motivated by more 
detailed models of galaxy formation such as semi-analytic 
models or hydrodynamic simulations (see e.g. Wechsler 
et al. 2001; Benson et al. 2001; White et al. 2001). Sim- 
ilar formalisms have been presented by e.g., Jing, Mo, & 
Borner (1998), Seljak (2001), Scoccimarro et al. (2001), 
and Berlind & Weinberg (2001). Simultaneous informa- 
tion on the number density of a population as well as its 
clustering on small and large scales can effectively con- 
strain all three parameters of the occupation function ex- 
pressed in this way (BWS02). However, measuring the 
small scale clustering requires accurate redshifts which are 
not available for the bulk of the EROs. Therefore, here we 
treat the slope of the occupation function as a free pa- 
rameter and investigate the allowed range of values for 
the minimum mass and normalization, based on the con- 
straints from the number density and large scale clustering 
of EROs. We use similar constraints for local giant ellip- 
ticals and the z ~ 3 Lyman Break Galaxies (LBGs) to 
derive their occupation functions, and use this informa- 
tion to speculate on the relationship of these populations 
to the intermediate redshift EROs. 

This empirical approach provides a useful complement 
to the detailed, physical models in several respects, per- 
haps particularly in the case of high redshift objects such 
as the EROs. Semi-analytic models based in the currently 
favored ACDM cosmology, which do well at reproducing 
the observed properties of z — and z — 3 galaxies 
(Somerville & Primack 1999; Somerville et al. 2001b; Cole 
et al. 2000), fail to reproduce the observed number den- 
sities of EROs, by a significant factor (Firth et al. 2002; 
Somcrvillc & Moustakas, in prep). The continuing uncer- 
tainty in the basic nature of the bulk of EROs (i.e. whether 
"old" or "dusty") makes it difficult to isolate the physics 
responsible for the discrepancy. The results presented here 
provide direct constraints on such models in terms of the 
basic theoretical quantities. Furthermore, in the hierar- 
chical paradigm, the age of the constituents of a system 
(e.g. stars) does not necessarily bear any relationship to 
the epoch of assembly of the object. Therefore, while nu- 
merous observations point to the conclusion that the stars 
in elliptical galaxies are rather uniformly old, this does not 
place any direct constraint on the mass assembly history 
of these galaxies. The clustering analysis that we present 
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here circumvents this ambiguity, again by directly con- 
straining the masses of the dark matter halos hosting the 
plausible progenitors of giant ellipticals at high redshift. 

The structure of the paper is as follows. In §2, we sum- 
marize the relevant observational data for the three pop- 
ulations that we will study: local giant ellipticals (gEs), 
intermediate redshift EROs, and high redshift LBGs. In 
§3, we provide a brief summary of the analytic model. In 
§4, we present constraints on the occupation function pa- 
rameters for the three populations, and in §5, we present 
predictions for the redshift evolution of the number den- 
sity and clustering properties of each population. In §6 we 
summarize our results and main conclusions. 

All computations are carried out for a ACDM cosmol- 
ogy, with f^A = 0.7, fio = 0.3, and a power-spectrum 
normalization of er 8 = 0.9. In calculating the power spec- 
trum, we assume H = 70kms _1 Mpc" 1 , but all results 
are scaled to hioo =H /(100kms _1 Mpc -1 ). 

2. OBSERVATIONAL CONSTRAINTS 

For any given population of interest, we require the co- 
moving number density and a measure of the clustering 
strength. In general, observed correlation functions are 
characterized by a power law functional form, 

<M-(£) T - (i) 

We can then convert this to a bias value, where the bias 
represents the clustering strength of the galaxies with 
respect to that of the dark matter (£dm) at a given epoch 
and scale. Several different definitions of bias appear in 
the literature; we use the square root of the ratio of the 
correlation functions at 8 hj" Q Mpc: 

b = (W£dm) 1/2 - (2) 
The correlation amplitude of the dark matter is straight- 
forward to obtain for an assumed cosmological model, us- 
ing cither numerical simulations or analytic approxima- 
tions. We have computed convenient fits for the dark 
matter correlation scale length ro and slope 7 from the 
GIF/ VIRGO N-body ACDM simulation (Jenkins et al. 
1998), which are described and shown in the Appendix. 
The dark matter correlation scale as a function of redshift 
is shown in Fig. 5. 

We consider three observed populations: 

• Local giant elliptical galaxies (gEs, z ~ 0) 

• Extremely Red Objects (EROs, z ~ 1.2) 

• Lyman Break Galaxies (LBGs, z ~ 3) 

The definition of each of these populations is somewhat ar- 
bitrary, particularly in the case of the high redshift objects, 
about which not very much is known. As we will rely on 
results from the literature, we are restricted to definitions 
used there. While the values of the number density and 
clustering statistics are inevitably sensitive to the details 
of how the sample is selected (especially to the magnitude 
limit), here we have chosen as a matter of philosophy not 
to address the uncertainties due to the ambiguity of sam- 
ple selection. Rather, we quote uncertainties that reflect 
the observational uncertainty associated with the determi- 
nation of the relevant quantity for a given sample. These 



values should be considered representative only. The re- 
sults should be updated as better observational data be- 
come available. 

Nearby "early type" galaxies may be defined in a num- 
ber of different ways, e.g. by morphological type, color, or 
spectral type. We wish to focus on luminous early type 
galaxies, as the EROs with which we want to compare 
them correspond to an intrinsically luminous population 
(if they are at the assumed redshifts z <~ 1). Locally, 
early type and red galaxies are well-known to be a factor 
of about 1.2-1.5 times more strongly clustered than the 
overall population at a given magnitude limit. Moreover, 
brighter galaxies also seem to be more strongly clustered 
than less luminous ones. However, as the bright end of 
the luminosity function becomes increasingly dominated 
by early type galaxies, it has not been clear up until now 
which effect dominates; i.e., whether the increased cluster- 
ing strength is due to the larger luminosities of early types 
or to the dominance of early types at bright luminosities. 
The recent results of Norberg et al. (2002) suggest the for- 
mer; i.e. that the luminosity dependence is the dominant 
effect. We therefore elect to use the values from Table 3 
of Zehavi et al. (2002), based on a volume-limited low- 
rcdshift sample of very luminous (—23.0 < M r * < —21.5) 
galaxies from the Sloan Digital Sky Survey (SDSS). While 
this sample has not been explicitly selected by type or 
color, it is more likely to be representative of the popu- 
lation in which we are interested than the fainter, color 
selected sample presented in the same paper. The values 
adopted (see Table 2) are similar to the results for L > L* 
early- type/red galaxies in the literature (e.g. Loveday et al. 
1995; Hermit et al. 1996; Willmer et al. 1998; Brown et al. 
2000; Cabanac et al. 2000). We expect that a more secure 
determination of the clustering of luminous red galaxies 
from the SDSS will be available shortly. 

The situation for EROs is at least as problematic. Un- 
fortunately there is no real consensus in the literature re- 
garding what should be deemed "extremely red," and in 
any case any simple color cut is bound to be somewhat 
arbitrary. This confusion is compounded by the variety 
of different filters used in these surveys. We refer the 
reader to Firth et al. (2002) for transformations between 
several commonly used sets of colors for typical galaxy 
spectra. Color cuts commonly used to define EROs in- 
clude R - K > 5 ot 6, R - H > A or 5, / - > 3 or 
4, / — K > 4. Another caution is that the surface densi- 
ties (and therefore number densities and spatial correlation 
scales) all appear to be very sensitive to the precise color 
cuts and magnitude limits assumed, and may be different 
for the "old" and "dusty" types. The values we adopt are 
drawn from the measurements of McCarthy et al. (2001) 
and Firth et al. (2002), which are based on galaxies with 
H < 20.5 and I — H > 3, corresponding approximately 
to a median rest-frame B-band magnitude of Mb = —20.3 
(Firth et al. 2002). 

Lyman Break galaxies (LBGs) arc defined as in the sam- 
ple of Adelberger et al. (1998), using a color-color cut and 
an (AB) magnitude limit of TZ < 25.5 (Steidel et al. 1996). 
We use the correlation function parameters determined 
by Adelberger (2000), and compute the comoving number 
density using the selection function given by Steidel et al. 
(1999); see BWS02 for details. We use the same overall 
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correction for incompleteness and the selection function 
adopted by Wechsler et al. (2001). 

The relevant clustering parameters for the dark matter, 
based on equations Al and A2, are presented in Table 1, 
and the observables discussed above are summarized in 
Table 2. 
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Fig. 1. — We plot the average bias and comoving number density 
as a function of minimum host halo mass M m i n , for the simple case 
of one galaxy per dark matter halo above the minimum mass M m i n . 
This corresponds to an occupation function with a = 0. 



3. AN ANALYTIC MODEL FOR GALAXY CLUSTERING 

Our model has only two moving parts. We adopt ana- 
lytic approximations connecting the masses of dark matter 
halos with their number densities and clustering strengths 
(Press & Schechter 1974; Mo & White 1996; Sheth & Tor- 
men 1999). We then parameterize the relationship be- 
tween dark matter halos and their galaxies (the occupation 
function) using a simple functional form. This allows us to 
predict the number density and large scale bias for a pop- 
ulation with a given occupation function, or conversely, to 
invert observed values of the number density and bias to 
obtain the corresponding occupation function parameters. 

For the halo mass function, we use the analytic expres- 
sion developed by Sheth & Tormen (1999), which agrees 
fairly well with the results of N-body simulations (see also 
Jenkins et al. 2001): 



dn h p da av 2 r 2 
dM = ^MuMYT" [ 



exp 



(3) 



Here, a is the linear rms variance of the power spec- 
trum on the mass scale M at redshift z and v = 8 c /a, 
where 6 C ~ 1.686 is the critical overdensity value for col- 
lapse. The other parameters are a — 0.707, p — 0.30, and 
c = 0.163, which were chosen to match N-body simula- 
tions with the same cosmology and power spectrum as the 
one we have assumed. The comoving number density of 
galaxies is then given by the integral over the halo mass 
function drih/dM, weighted by the appropriate galaxy oc- 



cupation function: 



dn h 
dM 



(M)N g (M)dM. 



(4) 



We determine the large-scale bias for galaxies by inte- 
grating the expected bias of halos as a function of mass 
bh{M), weighted by the galaxy occupation function N g : 



1 



L g ^M miI 
For the halo bias bh- 
Tormen (1999): 

b h {M) 



dn h 
dM 



(M)bh(M)N g (M)dM. 



(5) 



we use the expression of Sheth & 
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2p/S c 



(6) 



6 C l + {av 2 )P' 

The final missing piece is the galaxy occupation func- 
tion N g {M), or the number of observed galaxies per halo, 
at a given magnitude limit and redshift, and as a function 
of halo mass. In this paper, we choose a simple power law 
form, with a normalization Mi, a slope a and a low- mass 
cutoff M min : 



N g (M > M min ) = — 



M 
Mi 



(7) 



This functional form is a reasonably good approximation 
to the occupation function predicted by semi-analytic and 
hydrodynamic models (Wechsler et al. 2001; White, Hern- 
quist, & Springel 2001). Physically, the parameter M m i n 
represents the smallest mass halo that can ever host a 
galaxy, and Mi represents the mass of a halo that will 
host on average one galaxy. The slope a represents how 
strongly the number of galaxies per halo depends on halo 
mass. As one expects larger mass halos to host more galax- 
ies than smaller mass halos, we restrict our analysis to 
a > 0. The large-scale linear bias b g is independent of 
the normalization of the occupation function M\ and de- 
termines M m j n (for a fixed value of a; see Fig. 2). The 
normalization M\ is then fixed by the required number 
density. We also define the (galaxy-number weighted) "av- 
erage" host halo mass, (M): 

M^(M,z)N g (M,z)dM. (8) 



1 f°° 
(M) = - / 

n 9 J M min 



Once we have determined the occupation function pa- 
rameters using a set of observations at a fixed redshift, we 
can then calculate the observables (correlation length and 
number density) for populations at different redshifts, se- 
lected either assuming a constant minimum mass or con- 
stant number density (§4). This can be used to deter- 
mine whether observed populations at different redshifts 
are likely to occupy halos with similar masses — for ex- 
ample, to determine if the halos occupied by EROs are 
similar in mass to those that host nearby giant ellipticals. 
This is discussed in §5. 

4. SOLVING FOR THE GALAXY OCCUPATION FUNCTION 
PARAMETERS 

In this section, we use the observational constraints out- 
lined in §2 and the formalism summarized in §3 to obtain 
constraints on the allowed range of values for the occupa- 
tion function parameters M m j n and Mi. We carry out the 
analysis at redshift z = for giant ellipticals, at z = 1.2 for 
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EROs, and at z = 3 for LBGs. Subsequently, we will in- 
vestigate the implications for interpreting the connections 
between these populations. 

First, in Fig. 1 we show the average bias and number 
density for the simple case of one galaxy per halo, for ha- 
los above a minimum mass M m i n , as a function of M m j n 
(from Eqns. 4 and 5, with N g = 1). This corresponds to 
the case of a = in our model, in which case the parame- 
ter Mi becomes undefined, and the clustering and number 
density are solely determined by M m ; n . Note that for a 
given value of the bias, the maximum value of M m j n is 
obtained for this one-galaxy-per-halo (a = 0) case. From 
this, we may read off upper limits on M m ; n for 2 = gEs 

), 



(1.7 x lO^h^Mo), z = 1.2 EROs (1.5 x 10 13 hrnnM r . 



and z = 3 LBGs (7.7 x lO^h^Mg). These values are 
summarized in Table 3. It is interesting to note that the 
observed correlation length of present-day rich X-ray clus- 
ters is r ~ 20 hj"oQ Mpc (e.g. Collins et al. 2000), implying 
a very large bias of b ~ 10. From Fig. 1, we would infer 
masses of M m ; n £ few x lO 15 hj" o 1 o M for this population, 
consistent with X-ray and lensing masses. 
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Fig. 2. — The occupation function parameters for a fixed as- 
sumed slope, a = 0.8, at redshifts z = 0.0, 1.2, 3.0 (top, middle, and 
bottom panels). The constant-bias lines (vertical solid lines) are 
independent of Mi (which represents the mass of a halo that will 
host on average one galaxy). Combined with a measurement of the 
number density (curved dotted lines, in units of hf 00 Mpc -3 ) of a 
galaxy population, both Mi and M m i n are uniquely determined (for 
an assumed slope a). 
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Fig. 3. — The explicit inversion of the clustering and number density parameters for each galaxy population at its redshift recovers the 
minimum halo mass and the halo-mass weighted mass, M m i n (left panel) and (M) (right panel) respectively. This calculation was done for 
the full range of uncertainties in all observational quantities, as discussed in §2 and summarized in Table 2. Each point indicated corresponds 
to a different value for the slope of the occupation function, with the range a = — 1.0 in intervals of 0.1. A value of a f» 0.8 (shown by the 
double-circle points) is shown. The dot-dashed lines with triangles are for z = 0, L > L* Elliptical galaxies (gEs); the dashed lines with 
squares arc for z w 1.2 EROs; and the solid lines with diamonds are for z fa 3 LBGs. The labelled vertical lines mark the a = values for 
M m ; n (see Table 3). Note that the M m ; n and (M) results are relatively insensitive to the occupation function slope a for the strongly-clustered 
EROs and luminous ellipticals. The solid diagonal lines correspond to M m ; n = Mi and (M) = Mi. 



In Fig. 2, we show the solutions in the two-dimensional 
parameter space M m ; n vs. Mi for fixed values of the ob- 
served galaxy number density n g and bias b g , at the red- 
shifts of interest for our three populations. For purposes of 
illustration, we have fixed the slope of the occupation func- 
tion to a — 0.8. Whereas in the future it may be possible 
to use spatial clustering properties on small scales to con- 
strain the value of a (see BWS02), in this paper we must 
treat a effectively as a free parameter. We note, however, 
that in semi-analytic models (Somerville & Primack 1999; 
Kauffmann et al. 1999), the value of a is close to 0.7 — 0.9 
for bright galaxies, and so will concentrate on a fiducial 
value a = 0.8 in parts of the discussion. Fig. 2 allows us 
to read off the values of M m ; n and Mi for any desired values 
of the observables n g and b g . For example, at z = 0, using 
the values of n g and b g for L > L* ellipticals from Ta- 
ble 2, we can see these galaxies occupy halos with masses 
greater than M min ~ few x lO^h^M© (Table 3). This 
is consistent with the predictions of semi-analytic models 
(Somerville et al. 2001a; Benson et al. 2000), and with 
observationally determined group/cluster masses. 

In Fig. 3, we show the solutions for the occupation func- 
tion parameters for the three populations, for different val- 
ues of the slope a. Higher values of M m j n correspond to 
lower values of a. As noted above, the maximum value for 
M m j n is obtained for the limiting assumption of one object 
per halo, a = 0, and is depicted by labelled vertical lines 
for each case, for the central values of Table 2. The filled 
areas show the full range of parameter space allowed by 



the inversion of the chosen values for the number density 
and correlation length of each population with their stated 
(1 — a) uncertainties (Table 2). It should be noted that 
we have simply used the uncertainties from the literature 
for a given sample, and have not attempted to include the 
uncertainty /spread due to sample-to-sample variance or 
differences in sample selection. These effects will actually 
lead to a much larger region of allowed parameter space. 
The minimum mass M m j n is almost independent of a for 
the strongly clustered local ellipticals and EROs. For the 
more weakly clustered LBGs, the minimum halo mass is a 
stronger function of a. We also show the average host halo 
mass (M) (weighted by galaxy- number) , which is even less 
dependent on a for all cases. From this figure, we see that 
the clustering and number density of EROs are consistent 
with their occupying halos of similar masses to those that 
harbor present-day L > L* ellipticals. In contrast, LBGs 
must occupy halos that are several orders of magnitude 
smaller in mass. 

The relationship of the values of M m ; n and (M) to those 
of Mi is also interesting to consider. Recall that Mi cor- 
responds to the mass of a halo that harbors on average 
one galaxy — therefore halos with masses greater than 
Mi will tend to contain more than one galaxy. The case 
Mi » M m i„ (above the diagonal lines in Fig. 3) corre- 
sponds to a situation where most galaxies are the sole oc- 
cupants of their dark matter halos. The LBGs fall into this 
category. Conversely, for Mi ^ M m j n , most galaxies have 
one or more companions within a common dark matter 
halo. The EROs and gEs are in this category. This is inter- 
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esting, as it provides indirect evidence for the general con- 
cept of hierarchical clustering (i.e., that large structures 
like groups and clusters are built up over time through 
merging), and also is even more directly suggestive that 
EROs, like gEs, are found in group/cluster type environ- 
ments and may have been influenced by interactions with 
other galaxies. 

Although LBGs are not the main focus of this paper, as 
they have been the topic of so many analyses of this type, 
and because our results differ significantly from some of 
those in the literature, it is worth a brief digression to 
discuss the reasons for this difference. Some past analy- 
ses have used the clustering of LBGs to argue that they 
are harbored by halos with considerably larger masses of 
~ lO^h^MQ (e.g. Steidel et al. 1996; Giavalisco et al. 
1998). One reason for this is that early results on the ob- 
servational clustering values based on a few fields yielded 
considerably larger correlation lengths, tq = 6h^QQMpc 
(Adelberger et al. 1998), corresponding to bias values of 
b ~ 2.7 (using our definition of bias). Apparently, the 
LBGs in these fields showed an unusually strong degree of 
clustering, and more recent results based on considerably 
larger areas have yielded lower values; for example, Gi- 

S^iO^hrijMpc, 



± O.Sh^Mpc (see 



avalisco & Dickinson (2001) find ro 
and Adelberger (2000) finds tq = 3. 
the broader compilation and discussion in Wechsler et al. 
2001). The corresponding values of the linear bias (using 
our standard definition) are b ~ 1.4-1.9. Differing defi- 
nitions of bias lead to a still broader range of values; for 
example Adelberger et al. (1998) adopt b ~ 4 for the same 
ACDM model. In addition, most previous analyses of this 
kind have made the simplifying assumption that each halo 
above a given mass threshold M m i n contains exactly one 
galaxy, corresponding to the limiting a = case discussed 
above. As one can see from Fig. 1, had we assumed 6^4 
and a — for LBGs as was done in Adelberger et al. (1998) 
we would have found M min ~ 10 12 h 10 1 M Q , consistent with 
the value obtained by those authors. Note that our results 
are also consistent with those of Wechsler et al. (2001), 
in which a detailed treatment of LBG clustering based on 
semi-analytic modeling combined with N-body simulations 
was carried out, and with BWS02, which used an analytic 
model similar to the one used here. BWS02 obtained fur- 
ther constraints on the slope of the occupation function 



0.9 



1.1 from small scale clustering data 1 . 



5. RELATING POPULATIONS AT DIFFERENT REDSHIFTS 

It is always tempting to speculate on the connection of 
a population observed at high redshift with the more fa- 
miliar objects in the local Universe. This is a complicated 
and somewhat ambiguous task in the context of the hier- 
archical scenario, in which objects do not keep their iden- 
tity, but are constantly merging and can even change their 
Hubble type. We cannot try to "follow" the properties of a 
particular object through time with the simple models pre- 
sented here, but we can place some constraints on the con- 
nections of the three populations we have been discussing 
with one another, and we can ask specific questions about 



populations selected in certain ways at different redshifts. 
This is the topic of this section. We consider two types of 
models: a classical "galaxy conserving" model, and hier- 
archical merging models for objects selected with a fixed 
mass threshold, or with a fixed number density. 
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Fig. 4. — The redshift evolution of the comoving correlation scale 
ro (top); and the linear bias (lower panel), in the galaxy conserving 
model, for populations as indicated. The filled regions correspond 
to measured uncertainties, as described in the text and in Fig. 5. 
The dark solid line in the upper panel shows the correlation scale 
of the dark matter. 



5.1. The galaxy conserving model 

First, we consider a very simple null- hypothesis model 
for describing clustering evolution. If galaxies are assumed 
to form with some "bias at birth" and then to evolve with- 
out merging or changing their luminosity, their relative po- 
sitions will change solely as dictated by gravity. This is a 
fluid flow problem, and the behavior of galaxies as tracers 
of the underlying matter will obey the continuity and Eu- 
ler equations (Peebles 1980). Fry (1996) worked out this 
conceptually simple "galaxy conserving" model, wherein 
the linear bias evolves by 



b(z) = l + (b Q - l)/D(z). 



(9) 



Here, D is the growth factor of the universe for our adopted 
cosmology, and bo is the "bias at birth" of the popula- 
tion. For bo we can use the bias values at the redshift 



Although BWS02 use the same values as we do for the observed correlation function parameters of LBGs, they use a different definition of 
bias. They chose to assume that the slope of the dark matter correlation function is the same as that of the galaxy correlation function, rather 
than defining the bias at a fixed scale, as we have done. This is why the fiducial value of the bias for LBGs in BWS02 is 2.4 while we adopt 
bLBG = 1-9. 
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of observation, and then trace the evolution forward and 
backward in time again. Implicit in this model is that 
there is no merging or any loss of identity for any of the 
galaxy-particles (see also Tegmark & Peebles 1998). In this 
model, as the universe expands, the bias decreases mono- 
tonically and the correlation scale tends towards that of 
the underlying dark matter (Fry 1996). The results using 
this model are shown in Fig. 4. If LBGs evolved in this 
way, they would end up being slightly less clustered than 
z = gEs. The extrapolation of EROs' clustering prop- 
erties to z « suggest that they would correspond to a 
very strongly clustered local population, even more biased 
than L > L* ellipticals today. 

5.2. Hierarchical Merging models 

We now use the hierarchical framework presented ear- 
lier to produce predictions for the clustering properties of 
populations selected at different redshifts in two different 
ways. In the first case, we select objects with a fixed mass 
threshold (constant minimum mass, "CMM"). Here, 
we assume that all parameters of the occupation func- 
tion remain the same, and that halos above a fixed mass 
threshold (M m ;„) may host galaxies, according to Eqn. 7. 
We use the occupation function parameters obtained in 
§4 for each of the three populations, and run these val- 
ues forward and backward in time to predict clustering 
properties over a range of cosmological epochs. We also 
consider the properties of populations selected to have a 
fixed number density at all redshifts (constant number 
density, "CND"). Here, we assume that the slope and 
normalization of the occupation function remain constant 
over time, but that the minimum mass required to host 
a galaxy changes, such that the comoving number density 
of galaxies remains constant. This can be directly con- 
nected with observations. As discussed at the beginning 
of this section, we note again that neither of these mod- 
els is meant to represent the actual evolution in time of 
a particular object. Instead, we represent the evolution 
of statistically defined samples, which may be related to 
similarly defined observed samples. 

Fig. 5 shows how the number density or minimum mass, 
correlation length, and bias change with time for popula- 
tions selected in these two ways. We have kept the occupa- 
tion function slope fixed to a = 0.8, and we have assumed 
that the slope of the galaxy correlation function j g is con- 
stant with redshift. In the constant mass case, because 
halos of a fixed mass represent rarer fluctuations in the 
density field at earlier times, the bias decreases rapidly 
with time. In the constant number density case, the re- 
quired M m j n increases with time, so that the bias evolves 
much less rapidly, but still in the same sense. 

In the CMM case, the mass-selected ERO-like popula- 
tion run forward in time are predicted to be slightly more 
numerous and more clustered than the z = gEs. This 
apparently contradictory result is due to the fact that we 
found that for a = 0.8, EROs occupy slightly more massive 
halos than gEs, but Mi was considerably larger because 
of the higher number density we adopted for EROs. Con- 
current evolution in a, Mi, and/or "/ g could change the 
details of this prediction; and moreover the value of the 
number density is quite uncertain. Generally though, this 
result supports the idea that EROs could be the already 



mostly assembled progenitors of z = luminous giant el- 
lipticals. Conversely, the mass-selected LBG-like popula- 
tion is much more numerous and less clustered than either 
EROs or gEs when run forward in time to z = 0, because 
of the relatively low-mass halos they occupy in our model. 
These results suggest that if LBGs are indeed the progen- 
itors of gEs, they must undergo significant mass-growth 
through merging between z — 3 and the present. More 
speculatively, this perhaps suggests that early type galax- 
ies may have been mostly assembled by z ~ 1 , but were in 
much smaller pieces at z ~ 3 (or else their progenitors are 
not associated with the LBGs). 

In the CND case, the number-density selected ERO-like 
objects projected forward in time end up in much more 
massive halos, and therefore much more clustered, than 
z = gEs (but are much less massive and clustered than 
rich X-ray clusters). The masses and clustering properties 
of the number-density selected LBG-like objects projected 
to z = overlap those of the gEs. This may lead to the 
same conclusion as the above: that LBGs may evolve into 
gEs by z — 0, but must grow in mass considerably. 

6. SUMMARY AND CONCLUSIONS 

Extremely Red Objects are an intriguing population 
which may be the z~l progenitors of giant ellipticals. Re- 
cent results from Wide-Field Near-IR selected surveys with 
multi-band photometry have made it possible to obtain 
angular clustering and photometric redshift estimates for 
these objects, allowing the number density and real-space 
correlation length to be derived from the angular clus- 
tering and surface density (Daddi et al. 2001; McCarthy 
et al. 2001; Firth et al. 2002). These data provide im- 
portant clues about the nature of EROs by constraining 
the masses of halos that harbor them in a ACDM hier- 
archical model of structure formation. We have used an 
analytic model, combining a simple parameterization of 
the occupation function of dark matter halos with ana- 
lytic approximations for halo number density and bias, to 
constrain the masses of halos harboring EROs. We apply 
the same approach to z = ellipticals and z = 3 LBGs, 
and using this information, speculate on the connection 
between these populations and EROs. Our main conclu- 
sions may be summarized as follows: 

• EROs at z ~ 1.2 occupy halos with min- 
imum masses M m i n ~ 10 13 h^ 1 M Q , with a 
galaxy-number weighted average halo mass of 
(M) ~ 5 x IO^Ii^Mq — in other words, EROs 
are probably located in rich groups or clusters, 
at z « 1. These masses are lower limits for the 
subset of "old" EROs, if they are associated with 
larger correlation lengths. These results are rather 
insensitive to the assumed value of the slope of the 
occupation function a. 

• Elliptical galaxies at z <~ with L > L* occupy 
halos with very similar masses to those harboring 
EROs, lending support to the hypothesis that 
EROs may be the already almost fully assembled 
progenitors of present day giant ellipticals. 

• Recent results on the measured clustering of 
Lyman-break galaxies at z ~ 3 imply that they 
occupy much smaller mass halos than EROs or 



MOUSTAKAS & SOMERVILLE 



9 




01234 01234 
redshift redshift 



Fig. 5. — left: In the constant minimum mass model, the redshift evolution of the comoving number density (top panel); the comoving 
correlation scale ro (middle panel); and the linear bias (lower panel). Each filled region reflects the uncertainties in the values of the occupation 
function parameters (as shown graphically in the inversion figures of Fig. 3), propagated with redshift by the recipe presented in the text, 
right: In the constant number density model, the top panel shows the redshift evolution of M m i n ; the second and third panels are as in the 
left. The dark matter correlation scale (Eqn. Al) is shown as a solid dark line. These plots are computed for a = 0.8. 
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local ellipticals, with M min w 10 10 " 11 h 10 1 M Q and 
(M) w lO 11_12 h lo 1 o M0. This does not necessarily 
imply that LBGs are not the progenitors of 
present-day L > L* ellipticals, only that their 
halos must experience considerable mass growth 
via merging if indeed LBGs are to grow into giant 
ellipticals. 
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APPENDIX 

The correlation function of the dark matter in our chosen cosmology has been computed directly from the GIF/ VIRGO 
N-body ACDM simulation (Jenkins et al. 1998), and is shown in Fig. Al for several redshifts. Approximating the 
correlation function on linear scales with a power-law, we parametrize the comoving correlation length scale ro and the 
correlation function slope 7 as a function of redshift in a convenient analytic form: 

log 10 (hr o %m) - logi (5.00) - 0.00142 - 
0.419 log 10 (l + z)- 

l.lllog 10 (l + z) 2 , (Al) 
7o DM =-1.74 + 0.169*. (A2) 




Fig. Al. — left: The spatial correlation function, £dm(0> °f dark matter based on the GIF/VIRGO ACDM simulation, computed directly 
at redshifts 0.0, 0.5, 1.0, 2.0, and 3.0. The correlation strength decreases with increasing redshift, as shown in the legend. Linear regressions 
are computed for t^dm > 2-5h]~ 00 Mpc (to the right of the dotted line), and are shown as solid lines over the measurements. From these fits, 

relevant for the r = 8hyjQMpc scale considered in this paper, the correlation length (fo t DM) and power (7dm) are computed. Best-fitting 
polynomial fits are given in Eqns. Al and A2, and shown with the data in the right panel. 
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Table 1 

Dark Matter Clustering Properties 



Parameter 


z = 


z = 1.2 


z = 3 


^8,DM 


0.439 


0.173 


0.087 


r ll,DM Ohm M P C ) 


4.981 


2.513 


1.101 


7dm 


-1.736 


-1.517 


-1.229 



Note. — The correlation amplitude at 8h 10() Mpc (£s,dm)> correlation length (t-o,dm), and correlation function slope (7dm), for the dark 
matter, at the redshifts of interest, as calculated from Eqns Al and A2, based on the GIF ACDM simulation (Jenkins ct al. 1998). 







Table 2 






Observational Values used for Analysis 




Parameter 


gEs 


EROs 


LBGs 




z = 


z = 1.2 


z = 3.0 


n (h? 00 Mpc- J ) 


6.0 -10" 4 


(1.0±0.1)-10~ 3 


(4.7±1.0)-10" a 


r o,g (Wo M P C ) 


7.42 ±0.33 


9.5 ±0.5 


3.8 ±0.3 


7g 


-1.76 ±0.04 


-1.8 


-1.60 ±0.15 




0.88 ±0.07 


1.36 ±0.13 


0.30 ±0.08 


b 


1.41 ±0.06 


2.7 ±0.1 


1.9 ±0.2 



Note. — Number density (n), correlation length (?"o, g ), correlation function slope (7 S ), correlation amplitude at 8h 10 Q Mpc (§g l9 ) and bias 
(6) for the observed populations of galaxies discussed in the text. 



Table 3 

Characteristic Population Masses 



Population 


z 


a = 




a = 0.8 




M min (h^oMo) 


M min KoMq) 


M 1 (hr O oM ) 


(M) KMg) 


gEs 


0.0 


1.7- 10 ia 


2.7-10 i;i 


4.4 TO 13 


8.9 -10 13 


EROs 


1.2 


1.6 -10 13 


9.1 -10 12 


3.9 TO 12 


3.8 -10 13 


LBGs 


3.0 


7.0 -10 10 


1.3 TO 10 


6.0 TO 12 


5.5-10 11 



Note. — Characteristic mass values for each population discussed in the text. The M m ; n for the a = case (one galaxy per DM halo) is 
given, as well as the M m i n , Mi and (M) for the adopted a = 0.8 case. These latter values correspond to the points marked with double-circles 
in Fig. 3. 



